EE3310 Class notes 


Version: Fall 2002 

These class notes were originally based on the handwritten notes of Larry Overzet. It is expected that 
they will be modified (improved?) as time goes on. This version was typed up by Matthew Goeckner. 


Solid State Electronic Devices - EE3310 
Class notes 
Introduction 


Homework Set 1 

Streetman Chap 1 # 1,3,4,12, Chap. 2 # 2,5 Assigned 8/22/02 Due 8/29/02 

Q: Why study electronic devices? 

A: They are the backbone of modem technology 

1) Computers. 

2) Scientific instruments. 

3) Cars and airplanes (sensors and actuators). 

4) Homes (radios, ovens, clocks, clothes dryers, etc.). 

5) Public bathrooms (Auto-on sinks and toilets). 

Q: Why study the physical operation? 

A: This is an engineering class. You are studying so that you know how to make better devices and 
tools. If you do not understand how a tool works, you cannot make a better tool. (Technicians and 
electricians can make a tool work but they cannot significantly improve it. They, however, are not 
trained to understand the basic operation of the tool.) 

1) Design systems (Can you get something to work or not?). 

2) Make new - improved - devices. 

3) Be able to keep up with new devices. 

Q: What devices will we study? 

A: 

1) Bulk semiconductors to resistors. 

2) P-n junction diodes and Schottky diodes. 

3) Field Effect Transistors (FETs) - This is the primary logic transistor! 

4) Bipolar junction transistors - This is the primary ‘power’ transistor! 

By course end, the students should know: 

1) How these devices act. 

2) Why these devices act the way they do. 
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3) Finally, the students should gain a “manure” detector. This can be described as the ability to 
judge whether or not a device should act in a given manner, i.e., if someone describes a device 
and says that its operational characteristics are “such and such”, the student should be able to 
briefly look at the situation and say “maybe” or “unlikely”. (Only a detailed study can give 
“absolutely” or “absolutely not”.) 

Let us start the class by describing just what is a ‘semiconductor’. 

1) The conductivity of semiconductors occupy the area between conductors and insulators. This 
implies that the conductivity can range over many orders of magnitude. Further, the conductivity 
of semiconductors can be adjusted through a number of means, each related to the physical 
properties of semiconductors. Typical methods for adjusting the conductivity of a 
semiconductor are: 

a. Temperature 

b. Purity (Doping) 

c. Optical excitation 

d. Electrical excitation. 


2) Materials come from 
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3) In most semiconductor devices, the atoms are arranged in crystals. Again, this is because of the 
physical properties of the material. The structures of solid materials are described with three 
main categories. (This can and is further subdivided.) These categories are: 

a. Amorphous 

b. Poly crystalline 

c. Crystalline 

To understand the distinction between these solid material types, we must first understand the concept of 
order. Order can be described as the repetition of identical structures or identical placement of atoms. 

o 

An example of this would be an atom that has six nearby atoms, each 5 A away, arranged in a pattern as 
such. 
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If one where to pick any other atom in the material and find the same arrangement, then the material 
would be described as having order. This order can be either Short Range Order, SRO, or Long Range 
Order. Short-range order is typically on the order of 100 inter atom distances or less, while long range is 
over distance greater than 1000 inter atom distances, with a transitional region in between. 


We will now discuss each of the solid material types in turn. 

Amorphous solids are such that the atoms that make up the material have some local order, i.e. SRO, but 
there is no Long Range Order, LRO. (Materials with no SRO or LRO are liquids.) 

Crystalline solids are such that the atoms have both SRO and LRO. 

Polycrystalline solids are such that there are a large number of small crystals ‘pasted’ together to make 
the larger piece. 

For the purposes of this class, crystals, as we have said before, are the most important of these types of 
solids. Because of this we need to understand crystals in more detail. 

WE now need to introduce some basic concepts: 

1) The crystal structure is known as the LATTICE or LATTICE STRUCTURE. 

2) The locations of each of the atoms in the lattice are known as the LATTICE POINTS. 

3) A UNIT CELL is a volume-enclosing group of atoms that can be used to describe the lattice by 
repeated translations (no rotations!). This is further restricted such that the translations of the 
cells must fill all of the crystalline volume and cells may not overlap. In this way, the structure 
is uniquely defined. 

4) A PRIMITIVE CELL is the smallest possible unit cell. 

Often primitive cells are not easy to work with and thus we often use slightly larger unit cells to describe 
the crystal. There are four of very simple - basic - unit cells that are often seen in crystalline structures. 
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IT SHOULD BE UNDERSTOOD THAT THESE ARE NOT ALL OF THE POSSIBLE 
STRUCTURES. These structures are: 

1) Simple Cubic, SC 



2) 


3) 


Here a, b, and c are the BASIS VECTORS along the edges of the standard SC cell. 
Body Center Cubic, BCC 


Here the ‘new’ atom is at 
Face Center Cubic, FCC 
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Here the ‘new’ atoms are at (a/2 + b/2), (b/2 + c/2), (a/2 + c/2), (a + b/2 + c/2), (a/2 + b + c/2), 
(a/2 + b/2 + c). 

4) Diamond Lattice 

The diamond lattice is fairly difficult to draw. However, it is very important as it is the typical 
lattice found with Si, the leading material used in the semiconductor industry. 



A Diamond lattice starts with a FCC and then adds four additional INTERAL atoms at locations 
a/4 + b/4 + c/4, away from each of the atoms. 


Now that we have described a few of the simple crystal types, we need to figure out how to describe a 
location in the crystal. We could use our basis vectors, a, b and c, but it has been found that this is not 
the most advantageous description. For that we turn to MILLER INDICES. Miller Indices define both 
planes in the crystal and the direction normal to said plane. As we know, all planes are defined by three 
points. Thus, one can pick three Lattice points in the crystal and hence define a plane. From these three 
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points, we can find an origin that is such that travel from the origin to each lattice point is only along one 
basis vector and the distance is an integer multiple of that same basis vector. Thus our points are located 
at ia, jb and kc, where i, j and k are integers. 



To determine Miller Indices does the following: 

1) Determine the proper origin and the associated integers, i, j k. 

2) Invert i, j and k. Thus (i,j Jc) => ( l/i,l/j ,1/k). In our picture above we find that (3,5,3) goes to 
(1/3, 1/5, 1/3). 

3) Next, one finds the least common multiple of i, j and k and use this to multiple the fractions. In 
our picture above that multiple is 15. Thus we find (1/3, 1/5, 1/3) goes to (5,3,5). This is our 
Miller Index. 

If one of the integers is negative, it is denoted with a ‘bar’ over the number. Thus (-5,3,5) is written 
(5,3,5) 

Often, multiple planes are equivalent. These are denoted with curly brackets {}. In a SC lattice, all 
faces appear to be the same. Thus {1,0,0} O (1,0,0), (0,1,0), (0,0,1), ( 1 ,0,0), (0, 1 ,0), (0,0, 1). 

In addition to the planes, we might also be interested in a vector, i.e. moving in a given direction for a 
given distance. In this field, vectors are denoted with square brackets, [] . Thus a vector V= 1 .5a+lb => 
[1.5, 1,0] or more commonly [3,2,0], since we always want to move from lattice point to lattice point. 
Equivalent vectors are given with angle brackets, <>. 

Of interest is that the plane given by (x,y,z) has a normal of [x,y,z]. 

Example Problem: 

Q: What fraction of a SC Lattice can be filled by the atoms? 

A: Let us assume that the atoms are perfect hard spheres. This is an approximation known as the 
“HARD SPHERE” approximation. Further let us assume that the atoms are touching their nearest 
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neighbor. This is known as the “HARD PACK” approximation. Now each of the sides of the SC have a 
length of a. (‘a’ is not to be confused with the vector ‘a’.) Thus the volume of the cube is a 3 . Now we 
need to determine how much of each atom is inside the cubic volume. For this we need to look at our 
picture of the SC lattice. 



Let us look in more detail at the atom at the origin. 



We see that 1/8 of each atom is inside the cube. Thus the 

r 


8*l/8*volume of an atom = y nr 3 = 1 7t| 


— I = ina 3 . This 


the atoms is in ~ 0.52 = 52% . 


total volume of atoms in the cube is 
means that the fraction of the volume filled by 
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Chapter 2 Carrier Modeling 

Read Sections 9.1 and 9.2 

The late 1800s and the early 1900s set the stage for modern electronic devices. A number of 
experiments showed that classical mechanics was not a good model for processes on the very small 
scale. Among these experiments were the following: 

1) Light passed through two slits clearly shows an interference pattern. This means that light must 
be treated as a wave. However, light hitting a metal surface causes the ejection of an electron, 
which indicates a particle nature for light. Further, it was found that the energy of the ejected 
electrons depends only on the frequency of the incident light and not the amount of light. 

2) Electrons passed through two slits clearly show an interference pattern but they had clearly been 
found to be particles. 

3) In 1911, Rutherford established that atoms were made of ‘solid’ core of protons and neutron 
surrounded by a much larger shell of electrons. For example Hydrogen has a proton at the center 
with a electron orbiting it. However, classic electromagnetism combined with classical 
mechanics implies that the electron must continue to lose energy (through radiation of 
electromagnetic waves - light) and collapse to the center of the atom. Clearly this was not 
happening. 

4) A spectrum of radiation (light) is observed to come from heated objects that did not follow 
standard electromagnetism. [This radiation is known as ‘blackbody’ radiation.] A theory based 
on the wave nature of light was not able to account for this - in fact the theory predicted what 
was known as ultraviolet catastrophe - where by the amount of energy given off in the UV went 
to infinity. 

5) Hydrogen atoms (and all other atoms and molecules) were found to give off light at well-defined 
frequencies. Further these frequencies exhibited a interesting series of patterns that did not 
follow any known model of the nature of physical matter. 

6) Electrons shot through a magnetic field were observed to have an associated magnetic field. 
Further this field could be either ‘up’ of ‘down’ but no place in between. 

A rapid series of new models were developed which began to explain these observed phenomena. 

1) 1901 Planck assumed that processes occurred in steps, ‘Quanta’, and thus was able to accurately 
predict Blackbody radiation. 

2) 1905 - Einstein successfully explained the photoelectric effect using a particle nature for light. 

3) 1913 - Bohr explained the spectra of the Hydrogen atom by assuming a quantized nature for the 
orbit of electrons around atoms. 

4) 1922 - Compton showed that photons can be scattered off of electrons 

5) 1924 - Pauli showed that some ‘particles’ are such that they cannot occupy the same location at 
the same time (The Pauli exclusion principle). 

6) 1925 - deBroglie showed that matter such as electrons and atoms exhibited a wave-like property 
as well as the standard particle-like property. p = h/X = tik , where p is the momentum, h is 
constant (Planck’s Constant), X is the wavelength, k is the wavenumber 2k IX and Pi = h/2n . 

7) 1926 - Schrodinger came up with a wave -based version of Quantum Mechanics. 
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8) 1927 - Heisenberg showed that you could not know both the time and energy or the momentum 
and position perfectly at the same instant. Specifically 

ApAx > fi 

AEAl > h 

9) Etc. 

We will look at three of these in a little detail so that you the students have a little understanding of the 
principles involved. 

The Photoelectric effect 




It is found that the electrons emitted from can be stopped from reaching the collector plate by applying a 
bias to the collector plate. If one plots the bias required to stop all of the electrons, one finds a very 
simple curve 
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Einstein explained this by hypothesizing that light is made up of localized bundles of electromagnetic 
energy called photons. Each of these photons had the same amount of energy, namely hv, where v is the 

frequency of the light and h is a constant, the slope of the line, known as Planck’s constant. 

Sommerfield later proposed a model of a conductor that looks like 


hv E=eV appl ied 



<t> - work 
function 

Free electrons 
(Fermi Sea) 


Thus, one finds that the electrons in the metal are ‘stuck’ in a potential energy well. The photons then 
supply all of their energy to a single electron. The electron uses the first part of the energy to overcome 
the potential energy well, and the rest remains as kinetic energy. 

Bohr model of the Hydrogen atom 

Bohr’s model of the Hydrogen atom was perhaps the first ‘true’ quantum model. It does a wonderful 
job of predicting the then measured frequency of light emitted from an atom. (It misses some ‘splitting’ 
of the lines that later improvements to the experiments found and later improved versions of the model 
deal with correctly.) 
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The basis of the model is that the path integrated angular momentum of the electron, while in orbit 
around an atom, is in discrete states that vary as integer multiples of h. Namely, 
p e = mvr 

= nh/2n 
= nfi 


v = ■ 


nil 


mr 


We now have two other equations to work with 
The energy of the electron 
E=K.E.+ P.E. 


= \mv 2 - 


e 

Kr 


The centripetal force on the electron 


mv~ 


centripetal 


Kr 




Kmv 2 


22 


Ktm 


me 


From this we note that r is a function of n. For n = 1, ‘ground’ state, we find 


Kn 


me 


0.5 29 A 


where a 0 is the Bohr radius and is the smallest radius at which the electron orbits the proton in the 
Hydrogen atom. Finally plugging both velocity and radius into our energy equation we find the energy 
of the electron, 
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We see that the total energy of the electron is ‘quantized’ with the smaller quantum number having more 
energy. Again, we can look at the ground state, n=l, and find 

E X =R 

4 

_ me 
2K 2 h 2 
= -13.56eV 


where R is the Rydberg constant and is also the amount energy required to remove an electron from a 
Hydrogen atom. (This is a processes known as ionization.) This ionization potential ‘exactly’ matches 
the experimentally measured ionization energy. The energy emitted/gained between the states is 
‘exactly’ the energy of the photons emitted/adsorbed. (Better experiment showed that the model was not 
perfect but very close.) 


We can extend this model somewhat by assuming that the binding (electric) potential is due to all of the 
charges inside the outer shell. Then we get 

KftV 


Zme 2 


= 0.529 An 2 Z = 1 


E 


Bohr — 


1 Z 2 me 4 

2 K 2 ft 2 n 2 


= -13.56eV/n 2 Z=1 

where Z is the number of protons less the number of non-outer shell electrons. 


We can now graphically look at the energy and radius as a function of n 
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900 


0 



n 


If we look at the potential well the electron is trapped in, we see that the higher the energy, the higher 
the expected radius. 

In a true Hydrogen atom, the electron is trapped between the repulsive ‘strong force’ and the attractive 
electromagnetic force. The potential well that is created between these forces looks like 
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Energy (eV) 


The one major item that Bohr’s model missed is a splitting of the levels, or ‘shells’. This splitting is due 
to a splitting in the allowed angular momentum and particle spin (internal angular momentum) in each 
shell. Thus we find each shell given by a label n has an allowed set of angular momentums, given by a 
labels 1, and labels m, as well as spin given by label s. 

The overall requirements are 
n>l 
L<n-1 
-L<m<L 
s=±l/2 

The label T is often replaced with 1=0 => ‘s’, 1=2 => ‘p’, 1=3 => ‘d’, 1=4 => ‘f’, (and then follow the 
alphabet). Thus an electron in shell n=3, 1=3 can be labeled 3d. The higher the quantum numbers n and 
L, the higher the energy. This means that our picture of the potential well now looks like 



We can have up to 2(2L+1) electrons in that state because of the possible m’s and ‘s’s. We often add a 

superscript to our label to tell us how many electrons are in a given state thus 

3d => 3d 5 

or 3d 2 

etc. 


Usually the lowest energy states are the first to fill - 


This is in fact why the periodic table is the shape that it is. The Noble gases are on the right hand side 
and have completely filled - or closed - outer shells. The element on the farthest left will have a 
[noblejns 1 configuration, i.e. [He]2s‘ is Lithium while [Ne]3s* is Sodium (Na). 


At the close of the 1920, two versions of full fledge Quantum Mechanics were proposed, a wave version 
of QM by Schrodinger and a particle version, employing matrices, by Heisenberg. These are equivalent 
yet different and can be used to independently solve problems. For what little QM we do do, we will 
predominately use Schrodinger’ s version is this class. 
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where ^(r,?) is the ‘wave function’, a probability function for the particle/wave, 3, 


— , etc and 

3f 


V 2 = 3 2 + 3 2 + 3 2 . We note that the wave function implies that we can not know exactly when and where 
a ‘particle’ is located at. At best we get a general idea of were it might be. This is a very important 
concept as it leads to the idea of tunneling, which is very important in some modern devices. 


The unfortunate part of QM is that the equations are very hard to solve for any real physical system. 

The Hydrogen atom has been explored in detail this way but it would take us most of the semester to go 
through these calculations. As we are interested in understanding devices instead, we now look at some 
approximate models that will give us a feel for what is physically happening. 


The first model of a physical system that we will look at using Schrodinger’ equation will be a square 
well potential. We do this for two reasons, 1) it is a very simple mimic of the Hydrogen atom and 2) it 
is very similar to real devices that we can build. The potential is such that it is 



lO 0 <x<L 
[oo elsewhere 


We can now apply this to Schrodinger’ s equation 
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(What we have just done, is known as separation of variables. It is a standard method for solving a 
multi-dimensional differential equation.) Thus 
( 


“lA 2 +(v(r) “ £) f (r)=0 

Finally going to one dimension we find 
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.^l + {y{x)-E)y r) = 0 

Outside the well, the wave function must be zero, as the potential is infinite. (Or else the second 
derivative is infinite which is unphysical.) Thus we find 

£-a A 2 + i?j^(.r) = 0 0 <x<L 

i//(.r) = 0 elsewhere 

We start by looking at the 0 to L part and integrate twice to find 
^-d;y/(x)=-Ey/(x) 


\ » , , ±W2mEx/ft 

y/(x) = y/ 0 e 
or 

i//(.v) = A 0 cos (^JlmEx /fij + B 0 sin ^JlmEx /hj 

Now our wave function must be continuous in both zeroth and first order derivatives, so that at x = 0 we 
find A 0 = 0. (Remember y/(x) = 0 elsewhere.) Now at x= L y/(x) = 0 so that 

i//(L) = 0 = B 0 sin [JlmEL/Jij 


sl2mEL 

= Ml 

h 




n =0,1,2,... 


_ h 2 n 2 n 2 
2 mL 2 


r 

y/(x) = B 0 sin 

V 



Finally, we typically normalize the wavefunction to 1, so that our total probability is ‘1’. This is done 
by integrating 
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Prop = J \/ (x)y/(x)dx = 1 


li 



How is this related to our Hydrogen Atom? 

First the higher the value of n the higher we move up the sides of the potential well. Now, if we look at 
both positive and negative direction of our potential well around the core of the Hydrogen atom, we see 
a shape that looks like 



Where we have ignored the core area where the electron is not allowed. 

This, we can model Hydrogen in a way that is very similar to the above. Further, we would expect to 
see that higher energy states correspond to being higher up the potential well. Because of the shape of 
the well, we expect the more energetic electrons to orbit at a distance further from the core. This is 
indeed what we see. 


Day 3 

Homework set 2 

Chapter 3 # 4, 7 ,8, 9 Due Sept 10 th , 2002 


Recap 
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We have now learned a few of things 


1) We can know thing only with so much certainty. This is governed by the Heisenberg uncertainty 
principle 

2) We know that particles can act like waves and electromagnetic waves can act like particles. 
Further the wavelength/momentum relation is given by p=h/X. 

3) Quantum Mechanics does an excellent job of describing atoms as well as how individual atoms 
are structured. 

4) The Pauli exclusion principle states that two electrons can not occupy the same state in the same 
location at the same time. 


We will briefly look at that last item as is concerns our study of solid materials. There are few things 
that we need to note. 

1) Atoms have discrete energy levels caused by the potential wells around the nucleus. 

2) Solids are made up a large number of atoms. These atoms have energy levels as well, but the 
potential wells are adjusted by the fields from the nearby atoms. Here the Pauli Exclusion 
principle comes into play. 

This leads us to a new issue. We are dealing with atoms that are in close proximity to each other. What 
happens in such cases? Well let’s put two atoms close together and draw the total potential well. This is 
effectively what happens when two atoms are bonded together. 



Here we see that shells 3 and 4 above in each of the atoms ‘mix’ with the states in the other atoms. This 
would imply that if both atoms had state 3 filled, then we would have two identical electrons orbiting the 
two atoms. This cannot happen, rather we get a small splitting of the energy states. Because the 
potential is lower, the energy of one of the states is typically lower while the other state may be slightly 
higher. In general, the total energy of the combined state 3 is lower. We known this because if the 
energy was higher, the combined particles would try to go to a lower state, e.g. an unbound state. 

Further the average potential that the electrons are sitting in is lower. (This can also be shown with 
QM.) 

When there are a large number of atoms, say N atoms, we get an equally large number of splits in the 
energy band structure. Thus, it very common for a gas of a certain species to have a set of very well 
defined sharp spectral lines, while a solid of the same species will have very broad spectral lines. 
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We bring this idea up because we are dealing with solid-state devices. Thus the interaction of multiple 
atoms and atomic species is important to our understanding of this topic. How these atoms bond 
together is critical to the characteristics of the devices. 


We will now examine bonds between atoms. They fall into four main categories. 

1) Ionic NaCl and all other salts 

2) Metallic Al, Na, Ag, Au, Fe, etc 

3) Covalent Si, Ge, C, etc 

4) Mixed GaAs, A1P, etc. 


Ionic: 

The first of these types of material is related to the complete transfer of an electron from one atom to 
another. Cl for example would like to have a closed top shell and thus it takes an electron from the Na 
to produce a [Ar] electron cloud. Sodium on the other hand would like to give up an electron, so to also 
have a closed shell, in this case [Ne]. [THESE OUTER SHELL ELECTRONS ARE KNOW AS 
VALANCE ELECTRONS.] Both of these acceptor/donor processes provide lower energy states. This 
means that the two particles Na + and Cl are electrostatically pulled together or bonded. The electrons in 
question, are not shared by the atoms. Picture wise, this looks like 


e 



Na = [Ne] 3s 1 => Na + = [Ne] Cl = [Ne] 3s 2 3p 5 => Cl = [Ne] 3s 2 3p 6 = [Ar] 
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Metallic: 

The second of these comes in two forms. The first form has only a few valance electrons in the outer 
orbital. These outer valance electrons thus tend to be weakly bound to the atoms and are ‘free’ to move 
around. An example of this type would be Sodium, Na = [Ne]3s‘. 



WE WILL DISCUSS THE SECOND TYPE OF METALS BELOW. 

Covalent: 

In the covalent bond, two atoms share one or more valance electrons. In this way, each atom thinks that 
it has a closed outer shell. Because the outer shell is closed, these materials are typically insulators - 
although some might also be semiconductors. (This in part depends on the size of the atoms. The 
smaller it is, the more likely it is to be an insulator.) An example of this is Carbon, C=[He]2s z 2p 2 . 

e" 

e" c 

e 

e 
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C C C C c 



Mixed States: 

Mixed states are a combination of Covalent and ionic. An example is GaAs. Picture wise, these look 
like 


e 

e 


e — Ga 

e e 

Ga = [Ar] 4s 2 3d 10 4p‘ As = [Ar] 4s 2 3d 10 4p 3 

We see that Ga, which is column III, has three electrons in the outer shell while As, which is column V, 
has 5. As such the pair has 8 outer shell electrons, just enough to create a closed outer shell. Ga, it turns 
out wants to attract an additional electron more than As wants an additional electron. Thus one of the 
electrons spends more time near the Ga atom, making it partially negatively charged and the As partially 
positively charged. (A full electron is not transferred.) Thus, GaAs has some properties of covalent 
bonding and some properties of ionic bonding. The final structure looks like: 
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Energy Bands 

We can look at this in a second fashion. The properties of materials are the out growth of the splitting of 
the states in atoms that are close together. If we where to take N atoms and equally space them apart 
then slowly move them together, we would find that the splitting of the states grows as we get closer 
together. Hence for a single state we might see 


Energy 

level 



separation 

As our atoms get closer and closer, the more of the energy levels begin to split. Thus for an atomic 
species such as Si, ls 2 2s 2 2p 6 3s 2 3p 2 or [Ne] 3s 2 3p 2 , we get 
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3p 2 

2N electrons 
6N states 


3s 2 

2N electrons 
2N states 


If we could vary the separation of the Si, with in this diagram we see regions that correspond to two 
types of metal, a semiconductor and an insulator. (Note Si has a specific separation and hence it is a 
semiconductor.) 


Metal type l(seen at very large separations) 



Metal type 2 (seen at large separations) 



Semiconductor (seen at moderate separations) 
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Empty 


Narrow 

Forbidden 


Filled 


Dielectrics (insulators) (Seen at small separations) 
Allowed 

Empty 
Filled 



The allowed bands have specific names. The lower band is known as the Valance Band, as this is where 
all of the outer shell electrons will typically move, while the upper band is known as the conduction 
band as this is where we find conduction of electrons in solids. [NOTE EFECTRONS ARE NOT THE 
ONLY CHARGE CARRIERS IN SOLIDS. We will discuss this soon.] 


End here day 3 


What is important for conduction to occur, is that electrons must have ready access to allowed energy 
states that are empty. This is in the conduction band. This is because, for the electron to move 
physically, it needs to have both a position and an energy state to move into. By ready access, we mean 
that the electron must have enough available energy, through light or random motion (Think 
Temperature!) to be able to make the transition. What sort of energy might we be talking about? Well, 
room temperature is about 1/40 eV (or 1 eV ~ 11,000 K). Thus at room temperature, we might expect a 
significant number of the electrons to be able to gain ~l/40 eV. 

Aside on Temperature: 

The concept of temperature is relatively simple concept. If a material has a temperature that is 
above 0 K then there is some random internal motion. (This is very different than directed 
motion where <v>^0.) Often we find that the a materials motion is such that <v>=0, (velocity 
has direction and magnitude) while <lvl>?D (speed has only magnitude.) Further, because of the 
random statistical nature of atoms, the distribution of velocities is a ‘Normal’ or ‘bell curve’ 
distribution. (This is also known as a Maxwellian or Boltzman distribution.) The temperature is 
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a measure of the width of the distribution. Thus, the higher the temperature, the higher the 
variation in particle velocities. 

The normal distribution is 


p(y) = const —exp 
(7 


- {y-n ) 2 

_2 


Here a 2 is the population variance, <7 is the population standard deviation, and 77 is the central 
value. The distribution looks like 



y 


The Maxwellian distribution in terms of velocity (in 3-D) is given by 
f(v) = n^ 

While in terms of energy, it is given by 


r \3/2 

-m(v) 2 


\ 27tkT J CXP 

2kT 


f(<f) = n— exp 
kT F 


kT 


We will see this again soon. 


Now back to how Temperature influences conduction. We will look at C. If the material has enough 
internal (random) energy, some of the electrons in the covalent bonds may gain enough energy to break 
free. Picture wise this looks like: 
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Free 

electron 


Broken 

bond 


In terms of the energy band, it looks like: 


Energy 



Now we can go back and look at the concept of temperature in terms of the fraction of electrons that can 
jump from the valance band to the conduction band. 

We will use a Maxwellian distribution to approximate the number of electrons that might cross the band 


gap, 


kT 

? 

we see that the higher the temperature, the greater the chance that an electron will have enough energy. 


f(<f) = n— exp 
kT F 
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At room temperature, kT=l/40 eV. The band gap for carbon (diamond) is on the order of 3 or 4 eV. If 
we assume 4 for simplicity, we find that the fraction of electrons that are in the conduction band is 

Fract = J7 ^d^ 
n 


1 roo 

= kft g ex F 

) 

kT 

d^ 

= — kTexp 
kT F 

kT 


= exp 

[-41 

kT 



= exp[-160] 

= 3.3E - 70 !!!!! 


o 

This means that if the lattice constant for diamond is 4 A, then the number density of atoms is 
Number/volume = 8/(4A) 3 = 1.25E23 Carbon/cm 3 . Each carbon atom has 4 electrons in the conduction 
band and thus, we might expect about 1.6E-46 electrons/cm 3 in the conduction band. 

UNDERSTAND THAT THE ABOVE IS A ROUGH APPROXIMATION AND WOULD ONLY 
TRULY APPLY TO FERMIONS THAT DO NOT ‘INTERACT’ (Particles that do not interact but do 
follow Pauli Exclusion Principle). WE WILL GET TO A MORE APPROPRIATE METHOD 
SHORTLY. 


Carrier types and Carrier Properties 

From the above picture, we see that we have moved an electron into the conduction band. This means 
that it can move around and thus conduct current. However, we have also produced a vacant spot in the 
valance band. This means that the rest of the electrons in the valance band can now move but just into 
the empty spot - or ‘hole’. But if a valance band electron fills that hole, a new hole must be created 
from somewhere else. This allows conduction of current through the movement of ‘holes’ in the 
valance band. 

Now, if we apply an electric field to our material we get a force applied to both the conduction electron 
and the hole. That force is 

F = —eE (electron) 

F = +<?E (hole) 

Pictorially this looks like 
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Energy 



If we were in free space, we would expect 
F n = —eE = m„a (electron) 

F p = +eE = m p a (hole) 

where n and p stand for negative and positive or electron and hole. 

This brings up a valid question, What is the mass of a hole?!?!? 

To answer that question, let us first consider the electron. First, we are in a system that is clearly not 
free space. As the electron is accelerated toward the positive bias, it will collide with atoms. This will 
cause it to slow down. This means that the acceleration we see in free space will not be the acceleration 
we see in the material. We can solve this by defining an effective mass, m*, such that our force 
equation is correct. Holes will also have an effective mass, as they also move via the movement of 
valance band electrons. Thus we get in the material 

F )( = -eE = m*a (electron) 

F ; , = +eE - m* a (hole) 

[YOU MUST ALWAYS USE THE EFFECTIVE MASS IN YOUR CALCULATIONS!] 

Some times m*>m e and sometimes m*<m e . It depends on the specific system. 


Holes are not real particles but they act like them. They are in effect the average movement of electrons 
in the valance band. Because they are easier to deal with than all of the valance electrons, we keep track 
of them in the valance band. In the conduction band, the number of empty spots is much larger than the 
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number of electrons. Thus, in the conduction band, we keep track of the electrons in the conduction 
band and holes in the valance band. 


Types of Semiconductors 

In the above discussion, we have created electrons in the conduction band and holes in the valance band 
by liberating an electron from its bound state in the valance band. This means that the number of holes 
and electrons, or p and n carriers, is identical. If we lose one, it is through recombination with a particle 
of the opposite type. Hence, to get rid of an electron in the conduction band we need to move it into the 
valance band and have it fill a hole. (This sort of process occurs all of the time. Likewise, we have new 
electron-hole pairs being created all of the time. In equilibrium, the number of electron-hole pairs 
created in a time period is equal to the number destroyed.) 

Semiconductors that naturally operate with equal numbers of electrons and holes are known as 
INTRINSIC or pure (less common) semiconductors. 

There is also a second class of semiconductors that has an unequal number of holes and electrons. 

These are known as EXTRINSIC or DOPED semiconductors. There are two basic types of 
semiconductor dopants, n-type and p-type. n-type give rise to excess negative (electrons in the 
conduction band) charge carriers, while p-type give rise to excess positive (holes in the valance band) 
charge carriers. We will discuss each of these in turn. 


n-type and p-type dopants 

n-type dopants, also known as DONORS, operate in the following manner: 

For Si technology, donors come from column V. THUS WE MOVE UP ONE COLUMN TO GET A 
DONOR. The three most prevalent are P {[Ne]3s 2 3p 3 }, As {[Ar]4s 2 3d 10 4p 3 } and Sb {[Kr]5s 2 4d 10 5p 3 }. 
We will look at the valance structure of a P atom. 


e 



Now we will now place this in our Si lattice structure 
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Because this electron is an extra in the covalent bonds of the Si crystal, it takes very little energy to 
remove it from the P atom and thus move in to the conduction band. [This can be seen from the Bohr 
model of the hydrogen atom.] In terms of the energy band, it looks like: 


Energy 


Conduction 

band 


E d 


Donor 

state 


E v 


Valance 

Band 


► 


Position 


p-type dopants, also known as ACCEPTORS, operate in the following manner: 
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For Si technology, donors come from column III. THUS WE MOVE DOWN ONE COLUMN TO GET 
A DONOR. The most prevalent one is B {[Hc]3s 2 3p 1 }. We will look at the valance structure of a B 
atom. 


e 



Now we will now place this in our Si lattice structure 


Si — Si — Si — Si 


Si — Si — Si — Si 



Si — Si 


Si — Si — Si — Si 


Because this hole is an extra in the covalent bonds of the Si crystal, it takes very little energy to add an 
electron to the B atom and thus remove in from the valance band. [This can be seen from the Bohr 
model of the hydrogen atom.] In terms of the energy band, it looks like: 

[ONE OTHER THING TO NOTE IS THAT THE B ATOM IS SMALLER THAN THE Si ATOMS. 
THIS PLAYS A ROLE IN HOW THE CRYSTAL IS STRUCTURED.] 
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Energy 


E c 


E v 


Position 

ACCEPTOR, DONORS and AMPHOTERIC atoms in III-V semiconductors. 

As with Si/Ge types of technology, to create an acceptor type state, we need to add a dopant that has one 
fewer valance electrons than the base material. For GaAs, this means we note that our lowest column is 
Ga of column III and thus we need to pick something from column II, e.g. Be or Mg. This would then 
sit in the Ga lattice site and act as a hole. To create a Donor, we would have to pick something with one 
more valance electron and thus from column VI, e.g. S or Se. This would sit in an As site and act as 
donor. Amphoteric dopants are those that can be either acceptors or donors, depending on which lattice 
site they occupy. For example Si dopant replacing Ga has an extra electron and thus is a donor. On the 
other hand, Si on an As site is missing an electron and thus is an acceptor. Species with this dual ability 
are known as Amphoteric dopants. 
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[Again note difference in size of atoms.] 
End day 4 


Carrier Densities 

To get to the concept of carrier density, we first have to understand two things 

1) How the particles interact and what the appropriate probability distributions. 

2) The distribution of energy states available for the charge carriers. 

We will start with item 1 . 

Energy distributions 

There are in general two types of particles. 

1) Bosons 

2) Fermions 

These particle types are separated by whether-or-not they follow the Pauli exclusion principle. Boson 
do not interact and hence do not follow Pauli. [Further all Bosons have integer angular momentum 
states (spins), e.g. 0, 1,2...] Photons are an example of a Boson. Since Bosons do not interact we can 
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have as many of them in a place-time as we like. All Bosons follow BOSE-EINSTEIN statistics. 

[Where you think of Bosons, think of BOZO THE CLOWN. Clowns are infamous for being able to 
pack dozens of themselves, many more than imaginable, in very small cars.] 

Fermions do interact and do follow the Pauli exclusion principle. Further Fermions all have half integer 
spins (1/2, 3/2 etc.). Electrons are an example of Fermions. [Since Fermions are ‘Firm,’ you cannot put 
two of them in the same place at the same time.] Fermions can be further broken into two categories, 

a) Those with over overlapping wave functions - and hence always interacting 

b) Those without or only briefly overlapping wave functions - and hence seldom or never 
interacting. 

Type (b) is what we typically associate with gas molecules. We know that these follow a ‘bell’ shaped 
or normal distribution. This is known as a Maxwell-Boltzmann distribution. 

i \-c 

f(^) = n — exp 

kT LkTj 


Type (a) is what we are dealing with in solid-state devices, electrons that are constantly interacting and 
that obey the Pauli exclusion principle. These particles follow a Fermi-Dirac function. The function is 
set by rules that are not always simple to follow all the way through - but it can be done - see the back 
of the book, Appendix III. [The exact derivation is not useful for the purposes of this class.] 


m=- 


l + exp 


(£-£ f ) 


We know what The Maxwell-Boltzmann distribution looks like from above. What does the Fermi-Dirac 
function look like? [Note that we are not calling it a distribution!] 


UTD EE3301 notes 


Page 34 of 79 


Last update 12:18 AM 10/13/02 




- 2-10 1 2 

E-Ef (eV) 


This is very different than the Maxwellian distribution and it is entirely due to the fact that the particle 
cannot occupy the same energy state at the time. Note also that we cannot normalize the distribution 
such that the integral over all energies is one. This implies that we are missing a piece of the puzzle, 
which we will get to shortly. However, we need to look at a few other things first. 


First: 


For relatively large energy shifts ^ - ^kT 

very like much the Maxwell-Boltzmann distribution. 


f(<£ - <£q > 3kT) = 


1 


1 + exp 


Qg-4>) 

kT 


we find that the Fermi-Dirac function becomes 


1 


exp 


kT 


= exp 


kT 


This means that at high relative energy separations, we can approximate our function with the 
Maxwellian, and get reasonably good results. [This is in fact what we did with our example of 
conduction electrons in diamond above.] 


Second: 

For our situation r -0 is known as the Fermi energy and is written . Thus we have 
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m=- 


l + exp 


(t-tp) 


kT 


£ 


F is very useful. This is true for a number of reasons. 
1) When we plug ^F into our distribution, we find 

f(A0 = — 


l + exp 
1 


(<f F -<f F ) 


kT 


1 


1 + 1 2 

This means that there are an equal number of holes and electrons at the Fermi energy. 


2) We also find that the function is ‘symmetric’ around the Fermi energy. By ‘symmetric,’ we 
mean that the number of holes at an energy ^x below the Fermi energy is equal to the number of 
electrons at that energy above the Fermi energy. This can be proved very easily 

f(4 + ^f) = — 


l + exp 


((^X +^f)-^f) 
kT 


1 


l + exp 


((4)) 

kT 


exp 

(-(4)f 

kT 

1 + exj 

I (-(<?*))] 

kT 


l + exp 

(-(4))] 



kT 


1 

l + exp 

(-(4» 

1 + exp 

'(-(4))' 

kT 

kT 


= 1 - 


1 + exp 


Ho) 


kT 




1 + exp 


((-^x + O - 0 


kT 
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We can now use this to look at where the Fermi energy should lie for an intrinsic semiconductor. By its 
very nature, the number of holes in the valance band must equal the number of electrons in the 
conduction band. For simplicity, we will assume that all of the holes and electrons are at their respective 

band edges, and C c. From the above discussion, we find that the Fermi energy must be equi-distant 
from the conduction and valance band edges. (The distribution is ‘symmetric’.) [THIS IS NOT QUITE 
CORRECT AS WE WILL SEE LATER.] Thus if we look at a energy band diagram we see: 


Energy 



Position 

We can further complicate the picture by adding the Fermi-Dirac function. However, in our graph of the 
function above, we have energy along the horizontal axis. In our energy diagrams, we have the energy 
along the vertical axis. Thus, we need to rotate the function to get: 


Energy 



0 1 Position 
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This is often not done as simply knowing the temperature and the Fermi energy uniquely defines the 
distribution over which electrons will spread. Thus often times one only sees the Fermi energy on the 
band energy diagram. 

Likewise, we tell something about the semiconductor by the position of the Fermi energy. 

1) If there are more holes than electrons (p>n) than the Fermi energy must be closer to the valance 
band. 

2) If there are more electrons than holes (n>p) than the Fermi energy must be closer to the 
conduction band. 

This is entirely due to the ‘symmetry’ of the function and can be seen by moving the function relative to 
the band structure. 


STATE DENSITY 

Up to this point we have ignored the fact that the Fermi-Dirac function cannot be normalized to a 
physically reasonable number. The reason that this is so, is because we have not included the fact that 
we have a set of discrete states that the electrons can occupy. The density of possible states is 
determined by quantum mechanical rules that are not simple to follow all the way through - but one can 
come up with a convincing argument as to why it has to be what it is - see the back of the book, 
Appendix IV. [The exact derivation is not useful for the purposes of this class.] This density is: 


9 F? 

N(£)d£ = r47tk 2 dk = —=■ 

(2ti) 3 7t 2 


f *\3/2 


m 


*2 

V h ) 


9 

N(£)d£ = r27tkdk = -^~wd<C 


(2 nY 


^ 1/2 d^ 3-D 


2 - D 


71 /T 


N(^)d^ = — r2dk = ~—^—(C~^ 2 d(C 1-D 

(27t) i nh 

Note that = g(<0 i n some books. These functions look like: 
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If we take into account that the energies under consideration are relative to the band edges, we find 
slightly different densities in the conduction and valance bands. (The difference is a simple sign flip 
the square-root of the energy term.) For three dimensions they are 


N c (i’)d^ = ^ r 


v/2 

2 


/ * \ 3/2 


71 


m 


r 


£-£ c d£ 


Conduction band 


N y (£)d£ = — y 
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^jfC v - flAc Valance band 
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Now to get the distribution of states that are fill (electrons) or empty (holes) we need to multiply the 
Fermi-Dirac function with the state distribution function. 

■f *\3/2 

m 


n(£)d£ = f(£)N c (<r)d<!7=- 



1 

V2f 

1 + exp 

kT 

7T [ 


n z 




Electrons in the Conduction band 


p(£)d£ = (1 - f(<0)N v (<0d<C = 


1 + exp 


O^f) 


kT 


-n/2 


f *\3/2 

m 


7C 


\* 2 J 


- cd c Holes in the Valance band 


J / 


We can look at these again in terms of our band structures: 


Energy 



If we move the Fermi Energy up or down we get very different results 
Up O more electrons - p-type dopant 



Down O more holes - n-type dopant 
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We note that the charge carriers tend to bunch around the band edge. 


To get the total number of electrons and holes, we simply integrate over the whole range of energies in 
each band. 

n 0 = J|Q^g, r U£)N c (£)d £ Electrons in the Conduction band 

Po = (l _ f(£))N v (£)d£ Holes in the Valance band 

Here n 0 and p 0 represent the numbers in thermal equilibrium. 

The above equations are difficult to work with - and in fact cannot be solved analytically. We can 
however use a simplifying assumption, that we have relatively large energy shifts ^ c = £ ~ ^0 - 3kT 
the equations become significantly easier to deal with. 


n= J“ n(£)d£ = j~ f(<C)N c (<C)d<C = J' 
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= N c exp 
Likewise 


kT 
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( *, 3 3/2 


p-2 


m p kT 

K 2xh 


= N v exp 


exp 

( 4 - 4 ) 

kT 


( 4 - 4 )' 


kT 


N c and N v are known as the effective density of conduction and valance band states. 


Energy 



A final way in which we can write these are in terms of the ‘intrinsic’ energy, E ; , and the ‘intrinsic’ 
density, n,. The intrinsic energy is the energy half way between the conduction and the valance band. 
(In reality, it is the Fermi energy for the intrinsic material and hence it only has to lie close to the mid 
energy.) The density is found from the hole/electron densities at that energy. 


nj = N v exp 


(4 ~ 4 ) 

kT 


=> n = n j exp 


(4 - 4 ) 

kT 


nj = N c exp 


(4 - 4 ) 

kT 



1 

1 

& 

l 

iA 

=> p = njexp 

kT 


or we can multiply the two forms together to get 


nf = N v N c exp 


( 4 - 4 ) 


kT 


exp 


(4 - 4 ) 

kT 


= N v N c exp 


~^G 

kT 


= np 


Additional ideas from Streetman and Banerjee 

We have been drawing pictures of the energy band structure, showing position on the horizontal axis 
and energy on the vertical axis. First we note that position is by its nature, 3-D. However, we also need 
to realize that the momentum of an electron is independent of the location. (Or at least is potentially 
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independent of the location.) This means that we can add three additional coordinates, one each for the 
momentums in the three different directions. In QM parlance, the total energy of the system is given by 


1 2 

E = -mv + V(r,v) 
but 

p = mv 


E = jr I- V(r,p) 

2m 

When this is used to operate on a wave function (our electron’s probability function) we get 

k 2 

E = — + V(r,k) 

2m 

where the vector k is known as the ‘eigenvalue’ of the ‘eigenstate’. [It is just a fixed vector quantity that 
depends on the state that the electron is in.] If there is no potential shift due to the momentum then we 
get 

k 2 

E — V(r) = - 
2m 

We can plot this and find that the energy-momentum curve acts like a parabola in energy-momentum 
space. In other words (pictures?): 

Direct Band Structure 
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Sometimes, we get a potential shift due to the momentum, causing a shift in our parabolic curve. 


Indirect Band Structure 



In general we can and often do have peaks and valleys in our E-k space plots. They can and often do 
become very ugly, including having/giving different m* for each of the valleys. Streetman and Banerjee 
show a picture along those lines, which I will not try to duplicate here. 

This however brings up an interesting issue. The electrons in the conduction band typically occupy the 
lower part of the band - or in reality any valley in the conduction band. One can understand this by 
think about what happens to a bunch of balls randomly drop on an area of small hills. From our 
understanding of classical mechanics, we would fully expect that the balls to be somehow distributed 
among the valleys, with the largest number in the deepest/widest/lowest valley. (A very deep but very 
skinny valley may not allow any balls to enter, while a very wide but very shallow valley may not trap 
many balls either - so balance must be reached.) 

Fortunately we can use our model derived above to arrive at the number of electrons in the different 
valleys in the conduction band. We will show how this is done by example. 


Example problem 

We start this making a few simple assumptions. 

1) All objects have some amount of kinetic energy. 

2) All objects have some amount of momentum. 

3) The effective mass in each of these valleys can and often is different. 
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We know from above that the density of states is given by: 

*V /2 

N c (<Od£ = ^f \ J<C-£ c d£ 

it h" 


Let us assume that we have two valleys. Number one due to the [1 ,0,0] planes and number two due to 
the [1,1,1] planes. Thus we have for valley one has 2 equivalent planes and thus: 

r- / *\3/2 

N c i(<£’)d<£’= X y R 


N c i(<Od<f = X ^C-C a ic 

planes tt ft J 

K Ui J 

while valley two has 8 equivalent planes and thus for valley two 

r- r * \3/2 


N c2 (^)d^= I 4 hr 

planes tt y Tl 


2 4-2 

n l « 




<C-<Cr7d<C 


Multiplying by the Fermi function and integrating (assuming non-degenerate) we find 
o S 2rcm* 1 kT 1 

n 01 = 22 72 exp VT 


n O i 27tm; 2 kT f 72 ("-(^2 - <?f) 

n 02 - 8 2 -3 exp — 


we can now ask at what temperature they will have the same number of electrons. We find this by 
setting n 01 = n 02 . A few algebra steps gives (provided I have done these steps correctly!): 

kT= 42-4' 

A, 4=1 

2 m nl 
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Start lecture 5 

Homework Set 3 Chapter 3: 10, 11, 12 

Note: A good web site to look at: http://jas2.eng.buffalo.edu/applets/ 

At this point we should probably stop at look at what we have done. 


First, we found that the electrons (holes) where distributed over a range of energies. How they are 
distributed is given by the Fermi Function: 


m=- 


l + exp 


(£-£ f ) 


Second, we found that the electrons (holes) can only occupy certain states, as determined by QM. The 
available state densities are given by: 


N c (<Od<f 



€>£ e 


Conduction band 







Valance band 


Third, we found that multiplying the functions gives us our electron (hole) energy distributions: 

* ni i n t n\*i s ni in 1 a /2 1 / n n i n . • .1 m 1 


n(£)d£ = f(£)N c (<?)d<?=- 


l + ex„(^l 


Electrons in the Conduction band 


p(<C)d£ = (l-m)N y (<C)d£= 1 Ur „ \ — VZV 

1 + expV^l) * U J 


V2| m*Y 

— r- — J£. r - CdC Holes in the Valance band 


We however, are only interested in the total number of electrons in the conduction band (holes in the 
valance band) and hence we integrate over all allowed states. These integrals are none trivial but 
provided the Fermi energy is not close (>3kT away) to either the conduction or valance band we can 
approximate the number densities in equilibrium as 

r-(4-^ F )l 

n 0 ~ N c exp kT — - 


Po ~ N v exp 
where 


Cv-^f) 
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N c = 2| 


/ * \3/2 

1 m n kT ' 


2nh z 


J 


N v = 2| 


or 


/ * \3/2 

' nipkT ' 


2fth 

( 


N v ,c = 2l 


* \ 3/2 

m p nkT 


2nK 


= 2^1 
v 27 in 2 ) 


j 

3/2 (* \ 3/2 

m p,n 


m P 


y 


2.510E19 cm 


-3 


( * A 3/2 
m p,n 

m P 


at 300K 


At this point we realized that we have defined things in very general terms. However, we can modify 
our semiconductor by adding dopants, etc. and so we came up with a new set of values that are just for 
the basic (intrinsic) material, E F => E ; , etc. 

Tf-4/ 


nj = N v exp 


nj = N c exp 



=> n 0 = njexp 


kT 



i 

x ^ 

~Ti 

1 

=>P0= njexp 

kT 


or we can multiply the two forms together to get 


n f = N v N c exp 
= N v N c exp| 


kT 

~^G 

kT 


exp 


-(4 -A) 


kT 


= n oPo 

eliminating the need to even know 4 and 4 . 

Are we missing anything? We need to know n 0 and p 0 to be able to do anything. 

We know how to get and 4 from the temperature, m P and m n. All three of these are measurable 

values. We have ways to measure 4j, which we will discuss latter in the class, while 4 and 4 have 
been eliminated from our equations. Thus the only things that we do not have are our Fermi and 
Intrinsic energies. How do we get them? 


To do this, we need to look at our system again. When we have an intrinsic material, we expect that at 0 
K all of the energy sites in the valance band will be filled and all of the sites in the conduction band will 
be empty. As the temperature is increased, some of the electrons will jump from the valance band to the 
conduction band, creating electron-hole pairs. This however means that the number of electrons is 
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always equal to the number of holes. This concept is known as charge neutrality. It is guided by more 
than just counting, it is also guided by the fact that any large separation of charges will lead to strong 
electric fields that tend to pull the charges back together again. Thus, we set 


n ° = P° 

Now plugging in our distribution functions, and noting for this system = we find 
n oi = N c exp 


1 

A? 

1 

i 


‘(A-i’i)l 

kT 

= p 0i = N v exp 

kT 


eliminating terms, we get: 


exp 
N c_ ‘ 

1 

A? 

i 

aT 

1 

L kT J 

N v 

exp 

AT 

i 

A? 

i 

kT 

= exp 

(a v +< ^c - 2 

kT 




_(A 

+ < ^c) 


N v ] 


2 

+ 2 

nJ 

_(4 

+ ^c) 

kT, 

_l In 

( *\ 

nip 


2 

i 111 

2 

* 

l m nj 


3/2 A 


This means that the intrinsic energy lies very near the mid gap energy, with a slight offset due to the 
effective mass ratio of the electrons and holes. This offset is typically very small. 


Now, all we need is the Fermi energy. We understand that the Fermi energy is set by energy at which 
we would expect to have the same number of electrons and holes. If we are dealing with an intrinsic 
material we have a way to get at that number. However, we are often dealing with a material that has 
been doped and hence has either excess holes (p-type) or excess electrons (n-type). We need to 
understand how these dopants affect the Fermi energy in order to understand how to calculate the Fermi 
energy. 


The distinction is that we have added either donors or acceptors. However, we should still have charge 
neutrality - the electric field is a powerful force! This means that we need to add up our positive 
charges and set them equal to our negative charges. 


Po + N D - n o + n a 

where and ^ A are the number of ionized donors and ionized acceptors respectfully. (Think about 
this for a minute. A donor that is ionized has given up an electron, which is now moving in the 
conduction band. An acceptor that is ionized has pulled an electron out of the valance band, leaving a 
hole to move in the conduction band. Neither ion is able to move and hence neither is a charge carrier.) 
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For donor and acceptor states, Fermi-Dirac statistics determine whether or not the state is filled. If we 
look at an energy diagram, we see: 

Energy 

E c 
E f 

E v 


Position 

We note 

1) that E d > E F => f(E D )»l O almost all donors are ionized positively 

2) that E a < E F => l-f(E D )»l <=> almost all donors are ionized negatively. 

Thus we assert that the number of donor ions and the number of acceptor ions are that same as the 
number of donors and the number of acceptors. WE ARE IN EFFECT ASSERTING THAT THE 
SYSTEM IS NON-DEGENERATE! 

This gives us two equations, 

Po +Nd = n o +n a 

n ? = n oPo 

Combining them gives 
2 

~ “ + N D =n 0 +N A 
n 0 

u 

0 = n 5 +h 0 (n a -Noj-nf 

or 

2 

— e + N a =Po +N d 
P o 

u 

0 = Po + Po( n d -N A)~ n f 

These are both quadratic equations that we can easily solve: 
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The two solutions with the negative signs do not make physical sense as this will give us negative 
densities for the charge carriers. Thus we find 

Po = ( Na 2 N °) + ^(N D -N A ) 2 +4nf = njexp 
no = ( Np 2 Na -* + | A /(N A -N D ) 2 +4nf = n^exp. 


kT 


We can now use these to get the Fermi energy is terms of measurable quantities. 


Finally, there are a few simplifications that are often made to the above equations: 

1) N a , N d = 0 => intrinsic n 0 =p 0 =n ; . 

2) N d - N a » n ; => Doped n-type. 

Then 

A /(N D -N A ) 2 +4nf « ^(Nd - N a ) 2 = |N D - N a | = N d - N A 

(N D -N A ) (n d -n a ) 
n 0 2 + 2 _JN D _jn A 

and 

(n a -n d ) (n d -n a ) 

P0 2 2 
Noting the p 0 n 0 is still equal to n ; 2 we still get 


which is small compared to n 0 . This says that the electrons are our majority carriers while the 
holes are our minority carriers. 

3) N A - N D » n ; => Doped p-type. 

Then almost copying the above math we get 
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= t/(n d -n a ) 2 =|n d -n a | = n a -n d 


n o 


(n d -n a ) ( (n a -n d ) 
2 2 


and 


4) 


Po 


(n a -n d ) (n a -n d ) 

+ r -N a -JN d 


2 7 2 
Noting the p 0 n 0 is still equal to n ; 2 we still get 

5? 

Po 


n o 


which is small compared to p 0 . This says that the holes are our majority carriers while the 
electrons are our minority carriers. 

N a -N d ~ n ; => we have to use the full equations. 


Example: 


Question 

We have Si at room temperature (30 °C or 300 K) doped with 10 16 cm 3 B atoms. What is: 

1) n, 

2) p 0 ,n 0 

3) E ; -E f 

4) E f -E v 

5) Draw the band structure 


Answer 

1) We can look up the intrinsic density from table 3.17 (page 89) in Streetman. It is 
n ; =1.5 X 10 10 cm' 3 . 

2) The doping level is: 

N a = 10 16 cm 3 . 

n d = o. 

N a -N d » n,. Thus we can use the approximation: 

Po = N A -N D = 10 16 cm' 3 . 

n 0 = n ; 2 /p 0 = 2.25 E 20/10 E 16 cm 3 = 2.25 E 4 cm' 3 . 

3) We can now get E ; - E F from 


Po = nj exp 


( 4 = - 4 ) 


kT 


or 


^-^ F =-kTln 



= -0.026eV In 


1.5E10 ^j 

1E16 J 


= 0.35 eV 


UTD EE3301 notes 


Page 52 of 79 


Last update 12:18 AM 10/13/02 



4) We can now get E F - E v from 

— (^p-^y) 


p 0 = N v exp| ^ 
£,-£, = kT Inf 


-i or 


P° ^Noting that 

/ * \ 3/2 


N v c = 2.510E19 cm 


-3 


p,n 


m P 


and m P 


=( 


m 3 h /2 + 2m 


3 / 2\ 2/3 _ 


hh ) 


= 0.81 m„ 


(see appendix III) 


5) 


[Note that this mass ratio comes from addition of the three p 0 s (one for each direction) that have 
the same exponential terms. This is equivalent to the ‘averaging’ we have done above.] we find 

N v = 1.02E19 cm -3 

We can now plug this into our above equation to get: 
tCp - <C Y = 0.180 eV 

Energy 
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Carrier action - motion, recombination and generation 


Carrier action is driven by a large number of processes. We can look at the in a general sense by 
looking at the densities of our charge carriers in a small volume. The first thing that we should ask is 
how does this density change in a given time? Hence, we want to know: 

— =? 
dt 


and &=? 
dt 


We can work with either of these equations to begin to understand the process. The first thing to realize 
is that the density can change either by production/destruction of the carriers or by carriers moving into 
or out of the volume. We can take this into account by expanding our derivative, 
dn dn dn dx dn dy dn dz 

dt dt dx dt dy dt dz dt 


dn dn dn dn 

dt dx dy y dz 


dn 

= — + vlVn 
dt 

This is known as the continuity equation. It describes conservation of particles. If particles are not 
conserved, and here that is possible - we can have electron-hole recombination etc. - then the right-hand 
side of the equation is just the production/destruction, 
dn dn 


dt dt 


+ vlVn = f, 


production/destruction 


The first term on the right-hand side describes the local gain/loss of the carriers, while the second term 
describes the flow into or out of a volume. Note that gain/loss is not the same as production/destruction. 
With this understanding, we can begin to look at how the charges move, as well as how they are 
produced/destroyed. (There is a similar equation describing conservation of energy.) We will start with 
current flow. 


Current flow 

At this point, we can calculate n 0 , p 0 , E ; , E F , etc but how does that relate to our real use of 
semiconducting materials, current flow? It is after all the flow of current, and the ability to turn that 
flow on and off that makes semiconductors useful. To answer that question, we need to know how 
electrons and holes move in the lattice. 


There are three useful terms which describe how our charge carriers move in the lattice. They are 
conductivity (a), mobility (pi) and Diffusion (D). We will now discuss where each of these comes from 
and their physical meaning. 

Without any external force, we expect our charge carriers to have a random spread of energy and to be 
randomly spread across the device. This randomness is attributed to a Brownian motion type of process. 
Brownian motion is the observed motion of dust particles in the air. If one where to sit in a room with a 
bright light, one can observe the motion of dust in the air. This motion is random in nature with dust 
particles moving at a variety of speeds in a variety of directions. An individual dust particle might move 
in one direction for a while then turn around and head in the opposite direction latter. (This is also 
known as a ‘random walk’ or a ‘drunken sailor random walk’.) Pictorially this looks like: 
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For dust particles in air, the collisions are between air molecules and the dust. For our charge carriers, 
the collisions are with impurities in the lattice and the lattice atoms. Which of these two scattering 
centers is more important depends on the temperature. 


Now let us apply an electric field to the device. This means that we now have a force on the charge 
carriers. 

F = qE 

assuming that the electric field is only in the x-direction, we find: 

F = qE = qE x sf 


dt 

_ d Px 

dt 

The total change in momentum is then simply the sum over all particles, giving 

dPtotal _ y' d Pi 
dt j_| dt 


= £ qEj- 


i=l 


(internal) 


— nqE avera g e or external 

Just looking at this equation would imply that the electrons/holes will accelerate ‘forever’. However, we 
don’t see this sort of phenomena - think of the current in a wire - rather we see a constant velocity. 
(What we have described should apply to a metal wire as well as our semiconductor.) Additionally, 
what we really see is an average over all velocities. This means that we have to have a way to slow 
down the electrons/holes. That way is through scattering with our atoms. In this, we can think about 
throwing a BB at wall - the wall does not notice but the BB stops or bounces in the opposite direction. 
We can describe this deceleration using the average time it takes an electron/hole to lose the average 
particle momentum. 

d (Ploss) _ (p) _ Ptotal 

dt r nr 

where the average is just the total momentum divided by the number of particles. 

At some point we will balance the acceleration due to the electric field with the deceleration due to 
collisions. Thus, 
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d (p total) _ d (Pg ain } d (PlQSs) _ Q 

dt dt dt 



rearranging this, and writing the momentum in terms of the average velocity, we get: 
(p)=rqE = m*(v) 


(v) 


n 


E 


m 


= pE 

where ji is our mobility of the electron or hole. Note that in the book, q = electron charge, -e. Here q is 
simply a place holder and requires the insertion of -e for electrons and +e for holes. To get the current 
flow in the correct direction, see below, the mobility is always given as a positive number. Thus we find 

,, _ ( v )_ t M 


M-p 


M = M 

e m ; 


This is all well and good but we still need to current flowing through our device. 
For a single species of charge carriers, this is simply: 

J = nq(v) 

rlql 

= nq|lE = nq— 
m 


= oE 

where a is the conductivity. If we have both charge carriers 
J = nq(v) 

= q( n oPn+PoPp) E 
= oE 


= E/p 

where p is the resistivity. 


At this point, we want to complete this analysis. Most notably, we need to get the total current and the 
total resistance. (Remember, we measure things in bulk not in differential areas.) 

The total current is 
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1= J Jda 

cross 

sectional 

area 


The resistance is 
L 

R = J dR(x) 

0 
L 


= J 


0 


dV(x) 

I(x) 

Edx 

J(x)A(x) 


( , A(x) 

Assuming constant values, which is typical, we get 
pL 
A 
L 


R 


oA. 

Pictorially, this looks like: 



Electrons are inserted at x = L and removed at x = 0. 
Holes are inserted at x = 0 and removed at x = L. 


So now we need to ask just what does the mobility mean? 

P = * 

m 
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We know what the charge and the mass are, but what is tau? 

Tau depends on the number of scattering centers and the temperature of the lattice. 

1) The more scattering centers (dopant or defects) the more often collisions will occur. Hence, the t 
is small and so is // . 

2) The higher the temperature => the more lattice movement => the more scattering => the smaller 
t => smaller pi . 

At low electric fields, i.e. low velocities, the collision rate is constant. Thus pi ~ constant. However, at 
higher velocities, the collision rate increases, until increasing the electric field does not significantly 
increase the velocity. A graph of V vs E looks like: 


10 7 


10 6 


10 5 


At higher electric fields, the collision start to heat the lattice, causing more collisions. At some point 
most of the energy starts to get dissipated as heating of the lattice and not as kinetic energy of the 
electrons. 

At this point, we have found that an electric field can cause current to flow. But what does this mean for 
the shape of our energy bands? 

The answer is relatively straight forward, it bends them. We can show this very simply. Our bands are 
based on the energy of the electrons. Namely: 

Energy = kE.+ PE. 

How do you show how much is k.E. and how much is P.E.? 

Adding E g energy to an electron at E v - the edge of the valance band - will put the electron in the 
conduction band but with no additional energy, i.e. k.E. Any additional energy goes to k.E. It also turns 
out that the P.E. must have a reference energy, which will call E ref . Thus, the potential energy for the 
electrons must be 
PE.= E C -E ref 
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and the k.E. must be 
kE.= E-E c 
Likewise holes have a 

PE.= E re f — E v anc j kE.= E v -E N 0 t e: increasing energy is in the opposite direction for electrons 
(up) and holes (down) in our energy diagrams. 

Now 

E = -VV = — 3 X W (in 1 - D) 
but 

E c~ E ref =-q V 

U 

V = -i(E c -E ret ) 

SO 

E = -VV 

U 

E x =-d x V (in 1 - D) 

— qd x (E c - E ref ) 

= qd x ( E c) 

or 

= q3 x ( E v) 

or 

= q3 x ( E i) 

noting that E v , E c and E ; are uniquely related. What does this do to our bands? 



Position Position 


At this point let us assume that our system is set up such that the electric field exists across the device 
but we have not connected the wires at the ends of the semiconductor. In essence, we have our 
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semiconductor sitting inside a capacitor. At this point, we can show that under these conditions, the 
Fermi energy does not change. (This assumes that there is no current flow. If there is a current flow 
then E F does change!) 

Assuming that there is no current flow, there is still transfer of charge across the ‘device’ - it is just that 
the net or total charge transfer has to be zero. Let us divide our sample into two parts, part A and part B . 



For an electron (or hole) to move from A to B, an open spot needs to exist in B at the same energy.. 
Thus the number of electrons jumping from A to B is going to be proportional to the number of 
electrons (holes) in A times the number open stops in B . 
rate from A to B °c N A f A (E)2N B (l-f B (E)) 

The rate from B to A has to be similar, 
rate from B to A °c N A (l-f A (E))£N B f B (E) 

Because we do not have a net current flow, these must be equal and thus 
rate from B to A = rate from A to B 

JJ 

N A f A (E) 2N B (1 - f B (E)) = N A (1 - f A (E)) ZNflfB (E) 

JJ 

f A (E) 1(1 - f B (E)) = (1 - f A (E)) Zf B (E) 

ii 

(f A (E) - f A (E)f B (E)) = (f B (E) - f A (E)f B (E)) 

ii 

f A (E) = f B (E) 

Jl + jE-EFAj/kTj- 1 = _|_ e (E-E FB )/kT j -1 


Efa - e fb 

Showing that at equilibrium the Fermi energy is a constant. 
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Because of this, we can use E F as the reference energy. Thus by applying a bias, we can make our 
material go from a p-type to an n-type material as we move across the ‘device’ (or vise-a- versa). 



Energy 


E c 

E f 

Ei 

E v 




Position 


► 


We now have bent bands and thus move electrons/holes but we have no flow ! Does this make any 
sense? 

To understand why this might be true, we need to look at the world again. Let us consider a bowl of 
water that is standing on a table. Now let us put a drop of food coloring into the center. (This is a fun 
game for small kids and those of us whom would like to be young still!) What happens to the color? It 
spreads out slowly and somewhat randomly. In effect we have all of the dye molecules undergoing 
Brownian motion - each following a random path away from the initial position. The average over all 
particles is given its own name - ‘Diffusion’. Like many random processes the ‘average’ is reasonably 
well behaved and fairly easy to model mathematically. [There is a whole science that studies random 
processes, known as Chaos Theory.] 

In general, we can divide any volume into very thin slices, each of width <L>, where <L> is the average 
collision length of our particles. Let us look at the transfer of particles between two volumes, V n and 
V n+1 . (This model applies not only to our electrons and holes but also to just about any particle transfer) 



When a particle undergoes a collision, it will likely be scattered into the next volume. (This is because 
the width of each volume is so small, the particle will not collide again until it has moved <L>, or the 
width of the cell - putting it into the next volume.) A particle in V n has an equal chance of moving into 
V n _ , as V n+] . We now need to determine the number of particles per unit area that move from V n to V n+1 
in a given time. This is known as the ‘Llux,’ V. Let us fix the position of all of the particles except 

those in volume V n . On average, all of the particles in V n will undergo a scattering collision in time x, 
where x is the average scattering collision time. The number of particles in V n is simply the density of 
particles times the volume. Lrom this we get: 


UTD EE3301 notes 


Page 61 of 79 


Last update 12:18 AM 10/13/02 





Total flow 


v n -*v n+1 


' (n n A n (L)) 


where A n is the cross sectional area of volume V n . To get the flux, we just need to divide by the area. 

i(%( L )) 


n 


v n -»v n+1 - 


Now let us release the particles in V n+1 . 

r v n ^v n+1 = r v n ^v n+1 -^V n+1 ^V n 


j( n n( L )) j( n n+l( L )) 

T T 


, = 1 ( n n ~ n n+l )( L ) 

2 T 

Now, we note that <L> is very small. Thus, in 3-D 
3n ( n n-n n+ i) 

dx (L) 

and thus ! 


1 (L) 2 3n 

2 T dx 


in 1 - D or 



in 3-D 


= -DVn 

where D is the diffusion coefficient. Note that D is a positive number and thus expect the particles to - 
on average - flow from the higher density to the lower density volumes. (This is true for any diffusion 
process.) 


We have these net fluxes for both electrons and holes. 
r n = -D n Vn 

r p = - D p V P 

Because of this, we have net currents due to diffusion. 

J n — — q n /" n — q n D n Vn 

Jp — ^p^p ~ ^~p = "^IpDp^P 
J diff = Jn Jp — — ^p^p^P 


Now let us go back to our situation in which we have an electric field dragging the electrons and holes 
across our ‘device’. For that we calculated a current of: 
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J E =nq(v) 

= q(n|I„ +p|Ip)E 


( 


=q 


*q n ^ 

n— L ^- i + p- 

n 


q P 


A 


* 

'p ; 


E 


(We got this in our last class period!) 

Looking at each species independently we get 

Jp = q P pqp E -qp D p v p 
Jn — qn n M-n E 
J total — Jp Jn 

Now in our system that we started out with the total current has to be zero. (We don’t have connections 
to the wires!) Thus, 

J total — Jp Jn — 0 

In fact, because of particle conservation, this equation must hold for each species independently. 
Jn=q n n f l n E + qn D n Vn = 0 

U 

np n E = -D n Vn 
-D n Vn 


E = 


nPn 


and 


Jp = qppq P E -qp D p v p = 0 

PlipE = D p Vp 




E = 


D„Vp 


PM^P 

Now what are Vn and ^P? 
n = n . e [-(Ei-E F )/lT] 

Vn = — !V e [^ E i -E F) ,kT ]vE i 
kT 1 

= — VEi 
kT 1 
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p = „ i el< E >- E F) ,kT l 


Vp = i e [ (E ^ EF)/kT lvE i 


= — VEj 
kT 1 


So, what is ^ E i . . . .???? (We have gone a long way to get here!) Well we already came up with this 
during our last class period. 

E = -VV 




E x =-d x V (in 1 - D) 
— qd x (E c - E ref ) 

= q3 x ( E c) 

or 

= q3 x ( E v) 

or 

= q3 x ( E i) 

so 

E = qV(Ei) 


SO! 


E = ^Vn 
n(i n 



-o n r 


n^nl 




kT 

D n = 

M-n 


q n 

and 
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E = 


PM-p 

^P. 

PM'd 


Vp 




li 



kT 

q P 


These are known as the Einstein relations 


At this point we are able to move our charge carriers around by electric fields, giving rise to average 
velocities and hence mobilities. The charge carriers also move around on their own, through the random 
scattering based process known as diffusion. Further, we know that the holes and electrons are the result 
of another random process, energy sharing through collisional processes. This leads us to an interesting 
question, can we understand the generation and loss of individual electrons/holes and if so, will that help 
our understanding of how electronic devices work? 

Homework Set 4 

Streetman: Chap 4 # 1,2, 3,4 ,5 .6, 7. Assigned 9/26/02 Due 10/3/02 

Electron-hole loss and gain - Recombination and Generation 

Loss of electron and holes is in some ways the simplest of the two processes that we want to look at. 

We know that an electron cannot ‘just evaporate’. Charge and internal spin are conserved quantities. 
Electrons have both (q = -1.6E-19 C, s = ±1/2). Likewise, except under very special circumstances, 
mass is also a conserved quantity. Electrons clearly have that as well (even if it is a very small 
9.1 IE-31 kg). This means that to get rid of an electron, we must put it someplace else. Well let us think 
for a minute. Our ‘electrons’ are just those electrons that have made it up into the conduction band and 
our ‘holes’ are just places in the valance band that do not have electrons - remember ‘holes’ are an 
artificial construct that makes the math a lot easier. Thus, the only obvious someplace else is for a 
electron to ‘recombine’ with a hole, i.e. an electron in the conduction band loses energy and drops into 
the valance band, filling a hole. Hence when we lose an electron, we also lose a hole. 

Gain of electrons and holes must also follow the same pattern as lose - if we gain and electron, we have 
to gain a hole. Again, this is due to the fact that electrons have charge, spin and mass. 

We find that the creation/destruction of electrons and holes happens at the same time. Because of this, 
we refer to them as an Electron-hole pairs, or EHP. We can look at this process graphically. 
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Side note: This creation/destruction that we are considering is only for a ‘device’ that has 
already been built. Doping of semiconductor material can and does independently 
change the relative concentrations of electrons and holes. Thus by adding an n-type 
dopant, we can add an electron without adding a hole. Once the ‘device’ is built, 
however, gain or loss of holes and electrons is only through the process we are describing 
here. 

We can now look at how fast we might lose EHPs. We know that for an electron to move to the valance 
band, a hole must exist there. Likewise, for a hole to move into the conduction band (same process as 
above, just looking at it the other way around) an electron must be in the conduction band. Thus the loss 
rate must be proportional to both the number of electrons and the number of holes. Hence 
r = a r n 0 p 0 
2 

= a r nj 

where r is the recombination rate and a r is the proportionality constant. Under thermal equilibrium 
conditions, the lose rate must be equal to the thermal generation rate, g therma i. Hence 
^thermal ~ 1 — tt r n ()P() 

2 

= a r nj 

This means that if we raise the temperature and hence raise g thermah we find that the intrinsic density also 
increases. This is simply what we have found with the Fermi function, so it is not a big surprise. 


At this point we need to look for those processes that might give rise to recombination or generation of 
EHPs. Some example include: 
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detector 

Direct - Band to Band 


R-G 

center 



Indirect - Requires a R-G Center to change momentum. (Si and Ge are like this.) 


Fast 

electron 
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Auger (This is pronounced O’jhay - a French name - not auger as in com auger) - or impact ionization 

The process that dominates depends on the conditions and the material type. To determine which is 
most important, we need to better understand the second and third of these three processes because they 
have the additional requirement that momentum must be transferred for the process to occur. 

In chapter 2, Streetman introduced the concept of direct and indirect band structures. 



Direct Band Structure Indirect Band Structure 

Here the electrons and hole lie right on top of the lines. (There are no other acceptable (k,E) states. 
Remember that the horizontal axis is momentum and not position.) 

Why are these important? Well it turns out that while photons can carry away energy, hv, they cannot 
carry momentum, hk/27l. Thus there are typically very few electrons that can recombine directly with 
holes by the simple emission of photons. (There still might be a few but they are of very limited 
number.) 

On the other hand, photon absorption can and does continue to occur. While the photon does not bring 
in momentum, it can cause an electron to jump from the valance band to the conduction band provided it 
has more energy then is necessary to overcome the direct energy band gap, E a . This new electron in the 
conduction band then will undergo collisions with the lattice, transferring momentum, to reach the 
minimum energy in the conduction band. (These collisions serve to heat up the lattice.) Under certain 
conditions, the electron can still make the jump, provided the photon energy is greater than E g . This 
occurs because the photon excited the electron into a very brief temporary virtual state. (Remember 
Heizenberg!) Provided the electron then rapidly (instantaneously!) undergoes a momentum transfer 
collision, it can be indirectly transferred to the conduction band. If there is no collision in the length of 
time of the virtual state then the electron decays back to the valance band giving up the adsorbed photon. 
Because this is a three-body problem, i.e. the photon, electron and scattering center must all be there at 
the same time, it occurs less frequently - but it still can and does occur. Thus we find: 
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no transfer 
limited transfer 
significant transfer 


hvi 


< E g 

> E g ,< E £ 
>E„ 


So, how much light is absorbed? 


Let us set up an experiment that looks like: 



A certain fraction of the photons will be absorbed in any given length dx. We can measure how many 
by ‘measuring’ the intensity of the light. Thus 

I cc # of photons 
dl = l(x + dx) - l(x) 

A # of photons = -# of photons absorbed 
= -(fraction of photons absorbed) X (I per photon) X (initial # of photons) 

U 

dl = -yldx 

we of course can integrate this equation to arrive at 

I = Io e ~ > * 

where I 0 is the initial intensity of the light. 

This light produces EHPs as it is absorbed, which means that we have now increased the EHP 
generation rate. Thus we should expect n and p to increase. 

We can find out how much using a simple ‘conservation law’. 
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dn dp 
dt dt 

= Gen rate - loss rate . 

We know from above that the thermal generation rate is equal to the recombination rate. 
§ thermal — r — <^r n P 

although now we are not in thermal equilibrium and thus 
n*n 0 

P*Po 

plugging this into the above equation, we find that 


dp 


dt 


a r rij + light gen rate - 

a r np 

thermal 

thermal 

gen rate 

loss rate 


At this point, let us turn off the light. The light generation rate goes immediately to zero but the density 
of holes and electrons do not immediately go to zero. Thus we find 
dn dp 
dt dt 

= a r n\ - a r n(t)p(t) 

So what are n (0 and P(0? Obviously we can divide it into a constant part that is due to equilibrium 
densities and a part that is due to our extraneous sources - such as the light shining on the ‘device’. 
Taking those assumptions, we find 
n(t) = n 0 + Z\n(t) 

p(t) = po + Ap(t) 

Further, as we are not adding dopant, so 
An(t) = Ap(t) 

These are known as the excess carrier concentrations (caused by our light). We can now plug these into 
our time -dependent equation. 
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dn _ d(n 0 + An(t)) _ dzin(t) 
dt dt dt 

= dp = d(p 0 + Ap(t)) = ddp(t) 
dt dt dt 

= a r nf - a r (n 0 + An(t))(p 0 + Z\p(t)) 

= a r n 2 - a r (n 0 po +Po^ n ( t ) + n o^p( t ) + z ^p( t )^ n ( t )) 

= a r nf - a r |n 2 + p 0 An(t) + n 0 Ap(t) + z\p(t)z\n(t)j 

= - a r (poAn(t) + n 0 Ap(t) + z\p(t)z\n(t)) 

= - a r (po^n(t) + n 0 An(t) + z\n(t)z\n(t)) 

= - a r ((po + n 0 )An + Z\n 2 j 


dztn 

dt 


= - a r ((po + n 0 )An + Z\n 2 j 


or 


- «r((PO+ n o)^P + ^P 2 ) 


dAp 

dt 

These equations are tough to solve analytically. However, it is often the case that we have a low level of 
injection (creation) of charges. By this we mean that the increase in charge carriers is small compared 
to the total number of charge carriers. 

Z\n,Z\p « n 0 +Pq 

or 

Z\n 2 ,Z\p 2 « An(n 0 +p 0 ),Z\p(n 0 +p 0 ) 

Using these, our equations collapse to 

dAn - «r((Po +n oM n ) 


dt 


and 


dAp 

dt 


- a r ((p 0 + n 0 )z\p) 


These are clearly solvable. 
An = Z\n(t = 0)e -t/T 

Ap = An{t = 0)e“ t/T 
where 


a r (po + n o) 
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Now like everything else, we break this down into a result for intrinsic, n-type and p-type materials. For 
intrinsic materials we get the above equations. For n-type and p-type, we get some simplifications and 
hopefully a better understanding of what is doing on. 

n-type => n 0 »p 0 . 

1 1 

T — = 

«r(PO+ n o) «r( n o) 

Now 

T=>T p 

To understand this label, we need to look at the total density of each type of charge carriers 
n(t) = n 0 + An(t) = n 0 

p(t) = Po + Ap(t) = Ap(t) is possible 

therefore tau has to apply to the p density and not really to the n density. 

We can do the same thing for p-type materials to arrive at 
- 1 _ 1 

«r(PO+ n o) «r(Po) 
where 

n(t) = n 0 + An(t) ~ Z\n(t) is possible 
p(t) = po + Ap(t) - po 


We can do the same thing with indirect recombination (or trapping) or other similar processes. 
(We may or may not get to it in this class but you should understand that other processes are 
possible.) For trapping we get 
1 


1 


where c n and c p are the ‘capture coefficients’ and N t is the trap density. Multiple such processes 
occurring at the same time will result in a tau that is the sum of all of the taus. 


At this point, we want to go back to our system that we had been shining light on. We will turn on the 
light and leave it on. After a little time, the system will come in to a new steady-state density condition. 
(I do not want to call it ‘equilibrium’ as that is used for n 0 and p 0 . We clearly will not be at those 
densities.) Let us go back to our basic equations, 
dn dp 


dt 


2 

a r nj + light gen rate - 

a r np 

thermal 

thermal 

gen rate 

loss rate 


When we are in steady state, the time derivatives become zero. Thus, 
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dt dt 

— OC r Ilj + g 0 pt _ ®r n P 

= «r n ? +gopt - «r( n 0 + ^ n )(P0 + A V>) 

As before, we will assume that 
An- Ap 
so that 

0 = a r nf +g opt - a r (n 0 + Z\n)(p 0 + An) 

= a r nf +g opt - ct r (n 0 Po + Anp 0 + n 0 Z\n + AnAn) 
= a r nf +g opt - oc y {n} + Anp () + n 0 An + AnAnj 
= gopf « r (^n(p 0 +n 0 ) + Z\n 2 ) 

U 

gopt =«r(^ n (PO+ n o) + ^ n2 ) 


Again, we assume low-level injection. (This assumption will hold for all but a very few devices.) That 
is 

An,Ap « n 0 +Pq 
or 

An 2 ,Ap 2 « An(n 0 +p 0 ),Ap(n 0 +p 0 ) 
hence 

= Ap = g opt a ^ po+nQ j = gopt* 

So what is g opt ? It is simply the number of photons that are absorbed per volume per time; hence the 
number of EHP created per volume per time. We calculated this earlier. 


Well we have now set up a situation in which one of our fundamental equations no longer ‘works’ for 
the conditions of our ‘device’. We had derived: 


Po = n i ex P 


( 4 = - 4 ) 

kT 


n 0 = n^exp 


( 4 - 4 ) 


kT 


but we no longer have densities that are at the equilibrium densities! Now, we have 
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p = p 0 + zip = n iexp 


= ? 


kT 


+ zip 


n = ng + An = njexp 


kT 


+ An 


= ? 

At this point, we can define new two ‘Fermi’ energies, F p and F n , to get the same form of equation as we 
had before. These are known as Quasi-Fermi energies . 


p = nj exp 

( F P 4 'i) 

kT 

n = nj exp 

kT 

Example: 



We have a sample of GaAs (which is a direct semiconductor and hence we do need to worry about 
traps). Let us dope the ‘device’ such that: 


N„ 


TO 14 cnr, 


10 6 cm 3 , 


T n ~ lO' 8 sec, g opt = 10 20 EHP cnr s' 1 . 


gap 


1.43 eV. 


1) Find n 0 , p 0 

2) Find An, A P 

3) Find F p and F n 

4) Minimum light energy (in W/cnr 3 ) to do this. 


Answer: 

1) N a » n ; => p 0 = N a = 10 14 cnT 3 . 

n 0 = n; 2 /p 0 = 10 2 cnr. 

2) WE know 
An = Ap 

gopt =«r(^ n (Po+ n o) + zin2 ) 
1 

a r (p 0 +n 0 ) 


1 


a v = 


t (po +n o) 
1 


=> 


10“ 8 (l0 14 +10“ 2 ) 


10 6 cm 3 / s 


We can now plug this into the second equation to get 
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gopt = «r (^n(Po + no ) + An 2 ) 

= 10“ 6 (z\n(l0 14 + 10“ 2 ) + An 2 ) 


= 10 


20 


plugging this into the solution to the quadratic equation we get 
An = Z\p = 10 cm' l~g 0 pt T because Z\n«p 0 +n 0 1 

3) Now to get the quasi-Fermi Energies 


n = njexp 


1 

1 

iA 
1 

p = njexp 

-( f p-4)" 

kT 

| kT 




( F n-A) = - kTln ( n i/ n ) 

= -kTln(nj/ (ng + An)) 
-0.357 eV 


( F p-A) = kTln ( n i/p) 

= kTln(nj/(p 0 + Ap)) 
- -0.477 eV 


4) We can get the minimum energy by assuming that we are only providing sufficient energy for 
the electrons to jump across the band gap. Hence, we need all of the photons to have the band 
energy. This implies 


Energy = 10 20 * 1 .43 ^ * 1 .6X10 19 — 

* EHP eV 


3 

cm's 


23 


Watts 


cnr 


To gauge how much light this is, an incandescent light is about 5% efficient. Thus a 100 W light 
bulb puts out about 5 W of light energy. This would mean that you would have to have 
approximately 5 100 W light bulbs of light deposited in each cm 3 - a lot of light! 

Continuity Equation 

Now we need to back to the continuity equation. Remember that it describes conservation of particles. 

If particles are not conserved, and here that is possible - we can have electron-hole recombination etc. - 
then the right-hand side of the equation is just the production/destruction, 
dn dn 


i production/destruction 


— = — + v2Vm 
dt dt 

dp _ dp _ 

"~jj7 — jjjT ' 2. Vp — Iproduction/destruction 

Let us rearrange the equation 

U 

— — v ZVn + f production/ destruction 

— — V 2Vp + f production/ destruction 
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Noting that the first term on the right-hand side is simply the spatial derivative of the current density, we 
get 


^ ^ ^ + fproduction/destruction 

^ n fproduction/destruction 
On 

1 


dp 

at 


V 2qvp + fproduction/destruction 


1 


V2J p +f, 

q P 


production/destruction 


We can also begin to fill in terms for creation/destruction of our charge carriers. First, we have optical 
excitation - hence creation. Then we have thermal excitation and recombination. Thus, we get 


dn 1 

dt q n 


dp 


optical 

excitation 


thermal 

excitation 


recombination 


V2J n + g opt + - « R n(t)p(t) 


optical 

excitation 


thermal 

excitation 


recombination 


— = VIJ P + g opt + a R ni - a R n(t)p(t) 

ot qp 

Now, however, we have shown that the last two terms can be combined to leave a simpler single term 

— = — V2J n +g opt +a R nf-a R (n 0 +Z\n)(p 0 +Z\p) 
ot q n 

= — VlJ n +g opt -Z\na R (n 0 +p 0 ) 

^n 

— ^ 2Jn §opt — 

q n p ^ 

P VlJ p +g 0 pt + a R n ? _a R ( n 0 + zin )(P0 + zi P) 

q p 


at 


= -— v 2Jp +g opt -z\pa R (n 0 +p 0 ) 

q p 

= _ q^ VIJp+Sopt_ ^ 

At this point, we can now put in what we have derived for our currents 
Jn — qn n Mm^ qn^n^ n 


Jp = q P pqp E -qp D p v p 

Plugging these in gives 
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^ §opt 

at q n T 

= n(i n V EE + D n V 2 n + g opt - — 
jjf = ~ ~ V l(q p n ^ p E - q p D p Vp) + g opt - — 

= -pi^pV EE + D p V 2 p + g opt - — 

finally, we need to look at the left-hand side of the equation and note that the density is made up of the 
equilibrium part and the injected part and further note that the equilibrium part is a constant; thus 
11 

^ = nn„V 5JJ + D n V 2 4„ + g opt - ^ 

= -PH P V SE + D p V 2 Ap + g„ p , - ^ 


At this point we can use this rather ugly equation to gain a little more information. This new 
information will come in handy later. The new information is based on the following thought. We now 
have an equation that links shift in position with lifetimes of our carriers. Can we determine how far one 
of these carriers might travel before they are destroyed through recombination? Let us look at the 
simplest case. Let us assume that we have a constant source of our carriers at a single point (perhaps 
where a wire is attached to the material) and that we are just allowing them to move through the material 
on their own. Hence we are in steady-state and we are not applying any electric field. Further, we do 
not have any source of new carriers in the material, hence g opt = 0. Thus 
11 


dAn ^ „ j, An 

— = 0 = -D n V% 

fUo- + DpV**-f 

We can solve these equations to get 
An = D n rV 2 An 

= Cje+I^ + c 2 e + ^ / ^“* 

= C 1 e“l r l /Ln +C 2 e + I r l /Ln 

Ap = -D„rV 2 Ap 


C 3 e 


-I'I'VV + C 4 e +|r|/ V 5 i' 

H r l /L p , r +l r l /L P 


= C 3 e 1 1 p +C 4 e 

where Irl is the distance from the point source to wherever we are looking. We can determine the 
constant by looking at the results at both Irl = 0 and Irl = ». 
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Z\n|| r =0 = An 

= Cj +C 2 
M i = 0 

i|rl=oo 

= C 1 0 + C 2 °° 

=> Cl = An, C 2 = 0 


4p|| r | =0 = A P 

= C3 +C4 
Api , = 0 

1 I f poo 

= C 3 0 + C 4 o= 

=> C 3 = Ap, C 4 = 0 
Thus, we have found that 
n(r) = n 0 + Z\n(r) 

= n 0 + An(r = 0)e -r/Ln 
p(r) = po + Ap(r) 

= p 0 + 4p(r = 0)e r/Lp 

This means that our excess concentration of carriers drops off exponentially as we move away from the 
source of the carriers. 

Now, we can do back and plug this into our equation for the current. This will let us know how fast the 
charge carriers move from one place to the next. 

Jn — 

= q n D n v(n 0 + Z\n(r = 0)e“ r/Ln ) 

= -q n T^DnMr = 0)e“ r/Ln 

= -qn7^ D n Zin (r) 

Jp — — QpDpVp 

= -q P D P v (po + M r = 0)e r/Lp j 

1 — r/T 

= q P ^— D P ^p( r =°) e p 

= q P 7~ D P zi p( r ) 
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The final thing that we want to look at before we move on to real devices is what happens to our Quasi- 
Fermi levels. Remember that we have defined them as a convenience in our equations. 


p = nj exp 


F P 


kT 


n = nj exp 


(F„ - A) 
kT 


Plugging these into our full equations for the currents, 
Jn = qn n Hn E + qn D n Vn 


= q n n Pn E + qn D n Vn iCXp 


(F„ - 4) 

kT 


qn n M-n E 


q n Dn 

kT 


njexp 


(Fn-i’i) 


kT 




qn n M-n E Ahi n (^ E n qn E ) 
— q n nVF n 


Jp q P PFp E q P D p Vp 

— /t n pVFp 
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